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ABSTRACT 

A precise electromagnetic measurement of the sky coordinates and redshift of a coa- 
lescing black hole binary holds the key for using its gravitational wave (GW) signal to 
constrain cosmological parameters and to test general relativity. Here we show that 
the merger of ^ 10^~'^ Mq black holes is generically followed by electromagnetic flares 
from tidally disrupted stars. The sudden recoil imparted to the merged black hole by 
GW emission promptly fills its loss cone and results in a tidal disruption rate of stars 
as high as ^ 0.1 yr^^. The prompt disruption of a star within a single galaxy over a 
short period provides a unique electromagnetic flag of a recent black hole coalescence 
event, and sequential disruptions could be used on their own to calibrate the expected 
rate of GW sources for pulsar timing arrays or the proposed Laser Interferometer 
Space Antenna. 
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1 INTRODUCTION 

Recently, general relativistic simulations demonstrated that 
the coalescence of a black hole binary is accompanied by 
the anisotropic emission of gravitational radiation, causing a 
typical recoil of hundreds o f km s~ ^ for the black h ole rem- 
nant l|Pre torius 2005; Bak er et al.l [2006: Campan elli et af] 
I2OO6I ). Binaries of supermassive black holes (SMBH) are a 
generic consequence of galaxy mergers (tColpi fc Dottill2009l : 
lEscala et al1l2005l : iMaver et al]|2007l : ICallegari et al.ll2009^ ■ 
and the resultant gravitational waves (GWs) are potentially 
detectable with the proposed Laser Interferometer Space 
Antenna (LISA) or existing Pulsar Timing Arrays (PTAs) 
such as NANOGravB 

LISA would be most sensitive t o binary mergers 
with a total mass Mbh ~ 1O'^"^M0 (|McWilliams et all 
I2OIOI). whereas th e PTA sensitivity peaks at ~ lO^M© 
i Sesana et ai]|2009l ). PTAs have a signific antly poorer local- 
izatio n ability, with a typical uncertaintv (|Sesana fc Vecchid 
I2OI0I ) of ~ 40 compared to < 1 square degrees for LISA. 
These positional errors limit the precise determination of 
the luminosity distance to merging binaries. An electro- 
magnetic (EM) counterpart would greatly reduce the po- 
sitional error to sub-arcsecond scales, and also determine 
the redshift of the source, which would enable its use as a 
"standard siren" (independent of the cosmic distance lad- 
der) for precision measurements of the dark energy equa- 



tion ofstate_^M^^ Hughe!j[2005l : lArun et al.|[2009l : ISchutj 
l2009l : lBode et al.ll2010t ). 



^ http://lisa.nasa.gov/ 
^ http://nanograv.org/ 



For these reasons, prompt EM signals are of primary 
importance for realizing the full potential of GW cosmol- 
ogy. The proposed prompt EM signals have so far all as- 
sumed the uncertain presence of a circumbinary accretion 
disk prior to coalescence. Dissipation of GW energy in the 
disk mi ght result in a weak E M transient shortly after the 
merger (|Kocsis fc Loebll2008l ). re-equilibration of the inner 
edge of the disk could create an X-ray bright ening on a 
timescale of lO-lO'^ yr (|Milosavlievic fc Phinney i 2005), and 
shocks produced by the GW-induced recoil m ight gener- 
ate EM reverberations after the recoil (|Lippai et al. 200^ ) 
which may take ^ lO'^ years to dissipate as enhanced in- 
frared luminosity jSchnittman fc Krolik|[20o3 ). It is not ob- 
vious whether these EM signals could be distinguished from 
the much more abundant sources of temporal variability 
in single SMBH quasars. Moreover, the luminosity of any 
circumbinary disk is expected to be significantly reduced 
by the cavity associated with the decoupling of the binary 
from the inner edge of the disk in the final stage of i nspira l 
(IMilosavlievic fc Phinnevl2005l : ISchnittman fc Krolikil2008l ) . 
The disk is not expected to refill the cavity and return to its 
full luminosity f or a time ^ 7(1-1- z)(Mbh/10'^M (;))^-^^ yr af- 
ter coalescence (|Milosavlievic fc Phinnevll2005l ). On longer 
timescales, the portion of the accretion disk that remains 
bound to the recoiled SMBH is expected to be detectable 
as a kinemati cally and eventuall y a spatially offset quasar 
fcoeb 2007; Bonning et al.ll2007l: [Shields fc Bonnind l2008l: 
iKomossa et al.l I2OO8I : IComerford et al.l l2009l : [Shields et all 
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|2009| ). although its hfetime is limited by th e supply of gas 
that can remain gravitationally bound to it (jBlecha fc Loebl 
l2008h . 

Here we show that the tidal disruption of stars pro- 
vides a prompt EM flag that does not depend on the prior 
existence of a gaseous disk in the vicinity of the merging 
binary, and can result from mergers of gas-poor galaxies. 
Recent observations have demonstrated that tidal disrup- 
tion events (T DEs) have a gen e ric lightcurve and emis- 



sion s pectrum jGezari et al.ll2006l . |2008| . |2009| : iDonlev et al] 
l2002l : (Esauei et al.ll2007l ') that are distinguishable from nor- 



mal quasar variability. Moreover, we find that SMBH recoil 
results in a sequence of TDEs over a timescale of decades 
and potentially years, with a rate that is --^ 4 orders of mag- 
nitu de higher than the typical TDE rate in normal galax- 
ies (jPonlev et al.ll2002l ). The existence of TDEs accompa- 
nying SMBH mergers has been studied in the past only for 
long time delays (~ lO'^- lO'^vr) before lllvanoy et al.|[2005l : 
IChen et al.ll2009l ) or after (|Komossa fc MerrittlbOOSl ) the bi- 
nary coalescence event. While previous studies have focused 
on mechanisms to slowly feed stars into an empty loss cone, 
here we show that GW recoil will instantaneously shift the 
loss cone to a non-empty region of phase space. 



2 PHYSICS OF THE LOSS CONE 

A star will be tidally disrupted by a SMBH of mass Mbh if 
it passes within the tidal distance. 



(1) 



where m* and r* are the mass and radius of the star and 
7? is a dimensionless constant of order unity (|Diener at al.l 
[l99^. In our discussion we adopt rj = 1 and assume an ap- 
proximate main sequ ence scaling law of r« oc m°'^ , which 
implies rt oc m, '^'^^ (iMagorrian fc Tremain3 Il999l ). Tidal 
disruption does not occur if rt is smaller than the event 
horizon radius r^, in which case the star is swallowed in- 
tact by the black hole. For non-spinning black holes and 
solar mass stars, TDEs are therefore possible for Mbh ;$ 
IQ^Mq. a significant black hole spin can allow for (angle- 
depen dcnt) TDEs by SMBHs with A'/bh < 7 x 10*Mq 
iBeloborodov et al.lll9'?3 '). During a TDE, half of the star's 
mass becomes unbound, while the rest flows on Keple- 
rian trajectories, until the associated gas strea ms return to 
pericenter and collisiona lly shock each other (|Reej 1 19881 : 
lEvans fc Kochanekl Il989h . These gas streams return at a 
charact eristic mass infall r ate M oc t~^^^ (roughly speak- 
mg; see iLodato et all (|2009l ') for a more thorough treatment) 
and form an accretion disk whose blackbody emission peaks 
in the UV or soft X-ray with luminosities comparable to a 
supernova (|Loeb fc Ulmeill 19971 : IStrubbe fc Quataert 20091 . 
Other sources of emission include line radiat ion from the 
unbound debris (|Kasen fc Ramirez-Ruij 120101 ) ■ and a pos- 
sible brief period of super-Eddington mass fallback. These 
features are useful for differentiating TDE flares from su- 
pernovae or quasar variabili ty, and some h ave already been 
applied to candidate events (jKomossall2004l ). 

In a spherical galaxy with a stationary SMBH, a star is 
tidally disrupted if its orbital angular momentum per unit 
mass falls below a critical value. 



where we have approximated the relevant orbit as nearly ra- 
dial. Such orbits define the so-called "loss cone" . The rate of 
TDEs for a stationary SMBH is set by relaxation processes. 
Inward of a certain galacto-centric radius the loss cone is 
empty, but outside of it there is a "pinhole" regime where 
the rate of scatter in a nd out of the loss cone is grea ter than 
the orbital frequency l|Magorrian fc Tremainelll999l ). 

Over the long orbital timescale of stars, the impulsive 
GW-induced recoil of the SMBH remnant from a binary 
coalescence event yields a nearly instant aneous change in 
the black hole velocity relative to the stars (jO'Learv fc Loea 
2009). Viewed from the rest frame of the black hole, there 
is a sudden shift in the velocity of all stars, yielding a new 
loss cone defined by 



j' = |x X (v - 



- Vk)|^ < Jcrit 



: 2GMBHrt, 



(3) 



where Vk is the SMBH kick velocity. 

We parametrize the density of stars around a SMBH 
binary in the last stages of its inspiral, using a power-law 
profile. 



p = poir/roY 



(4) 



This density profile corresponds to an isotropic pre-kick dis- 
tribution function of stars. 



2N7-3/2 



f{r,v) = C{2GMji^/T-v') 
where the normalizing constant is given by 

(3 - 7)(7 - 0.5)r(7 + 1) Mbh f ri^f 



C = 



(5) 



(6) 



27r2-5r(7 + 0.5) rf„f V2GA/bh 

which in turn depe nds on the observationally calibrated 
iMerritt et al.ll2009l ) radius of influence of the SMBH (in- 
side of which the mass in stars is 2A'/bh). For a full sample 
of galaxies this is 



rinf ^ 



24(Mbh/1O*M0)°-''' pc. 



while for core galaxies alone it is 



rinf 



■ 35(A/bh/1O^M0)°-'*'' pc. 



(7) 



(8) 



Stars will be bound to the black hole system of total mass 
A/bh before the binary coalescence if 



< 2GMw/r, 



and after coalescence if 



(v - Vfe)2 < 2GMw/r. 



(9) 



(10) 



J = |x X vr < Jcrit « 2GMBHrt, 



(2) 



The intersection of these two spheres in velocity space, with 
each other and with the loss cone, is the region of phase- 
space containing bound stars which are tidally disrupted 
after the recoil. By performing a Monte Carlo integration 
of the appropriate distribution function over this region we 
calculate the number of post-recoil TDEs. Through a sepa- 
rate integral we find that the unbound stars provide < 10% 
of the total number of TDEs and can be neglected. Another 
small correction to the total rate involves the net SMBH 
mass loss by G W emission, which is typ ically < 5% of the 
pre-merger mass (jCampanelli et al ] |200g ). We have included 
the associated small reduction (by ~ 10%) in the number of 
tidally disrupted stars. To evaluate the observability of the 
recoil-induced TDEs, we define the quantity N^{t) as the 
number of stars in the post-recoil loss cone which are tidally 
disrupted in t years. Since stars which fall into the new loss 
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Figure 1. Expected number of stars disrupted in less than 100 
years, A'^< (100), in the dry merger, joint core-cusp case. Left panel: 
Black hole mass dependence for kick velocity vj; = 400 km s~^. 
The dotted blue line (labeled 'BLC') represents the physically 
realistic scenario including the binary loss cone discussed in the 
text, whereas the black line (labeled 'NBLC') accounts only for 
the loss cone of a single black hole. Right panel: Kick velocity 
dependence for a black hole mass A/bh = IO^'^A/q, using the 
same line types. 



cone are near their apocenter in the rest frame of the kicked 
black hole, Af<(100) is the number of stars with orbital pe- 
riods below 200 years. 

The result of the Monte Carlo integral is sensitive to the 
innermost region of the distribution function, whose details 
are determined by the pre-merger dynamical environment. 
In particular, the SMBH binary will excavate a larger and 
more complex loss region than is given by Eq. ((3|. Below, 
we consider the innermost regions of phase space in gas-poor 
and gas-rich mergers, separately. 



3 DRY MERGERS 

In the absence of gas, 7 = 1.75 is the dynami- 
cally relaxed (so-calle d "Bahcall- Wolf " ) equilibrium state 
l|Bahcall fc Wolj|l976l ') of a stellar cluster around a SMBH. 
However, core galaxies are believed to be the end product 
of a binary inspiral, as the binary sheds angular momen- 
tum by ejecting stars. Numerical simulations of this pro- 
cess show that a binary hardening its orbit thro ugh scat- 
tering of stars will scour a core ijMerrittl l2006l l. but at 
some point depletion of the remaining stars in the binary 
loss cone (including in this context stars whose pericenters 
fall within twice the binary semimajor axis) will lead to a 
stallin g of the binary and the so-ca lled "final parsec prob- 
lem" jMilosavlievic fc Merrittll2003l ). Without gas, the bi- 
nary can only merge via a re-population of i ts loss cone. Sig- 
nifica nt triaxiality of the galaxy potential l|Merritt &: PoonI 
|2004 ) tends to re-populate the binary loss-cone but pre- 



serve a core with 7 ~ 1. Alternatively, collisional processes 
generate a central cusp of stars, though this method of bi- 
nary hardening is only effective for binary masses < 1O^M0 
(Me rritt et al. 2009,), and it may be modified in the presence 
of massive perturbers l|Perets fc Alexandeillioosh . such as in- 
falling molecular clouds in gas-rich mergers. These gas-free 
scenarios lead us to consider both core galaxies where 7=1, 
and galaxies with a joint core-cusp density profile where an 
inner 7 — 1.75 profile meets an ou ter 7 = 1 profile at a ra- 
dius of 0.2nnf (|Merritt et al.ll2009l ). We use Equation® as 
the radius of infiuence in these models. 

It is necessary to exclude the stars located in the pre- 
coalescence loss cone. The size of this region of phase space 
is somewhat uncertain. For collisi onal re-population of the 
loss cone, numerical simulations (iMerritt et al.ll2007^ indi- 
cate the SMBH binary will decouple from a relaxed distri- 
bution of stars at a semimajor axis of acq ~ lO'^rinf, and 
coalesce before the stars can relax into the gap left behind. 
To model this cavity in energy space we remove all stars with 
semimajor axes less than this radius. However, because re- 
laxation in angular momentum is faster than in energy, the 
resultant gap in angular momentum space will be partially 
refilled prior to the merger. The timescale f or filling up a gap 
in an gular momentum space is given by (|Merritt fc Wang 

[iooi), 



(11) 



where a is the semimajor axis of the SMB H binary, taken to 
be th e pericenter at which stars are ejected (| Merritt fc Wang 
|2005| ). a nd the system relax ation time at the radius of influ- 
ence is iMerritt et al.ll2009l '). 



T. ^ 8.0 X lOVr f 



(12) 



Thus, a second cavity in the distribution function is created 
by removing all stars with pre-kick pericenters lower than 
the binary separation a at which Tgap equals the gravita- 
tional wave timescale, 



Tgv 



256G3M|H^i' 



(13) 



with n — M1M2/MBH being the reduced mass of the binary 
and Mbu = {Mi + M2). For simplicity, we adopt a flat J 
dependence for f{E, J) with the cuts mentioned above. In 
the classical loss cone calculation, the steady state solution 
of the orbit-averaged Fokker-Planck equation yields a distri- 
bu tion function that varie s logarithmically with J at fixed 
E (|Cohn fc Kulsrudlll978l ). However, this solution may not 
apply to the innermost stars, for which the orbit averaged 
assumption may break down and strong star-star scatter- 
ings (which the Fokker-Planck approach does not account 
for) could be important. Furthermore, the loss cone of a bi- 
nary SMBH is not the pure sink assumed for a single BH, 
since stars may remain bound to the binary on low- J orbits. 
The use of a logarithmic instead of a flat distribution would 
have reduced A'^< (t) by a factor of ~ 2 — 4. 

We use these cuts in stellar energy and angular mo- 
mentum as modifications to the joint core-cusp profile. We 
generalize the results to a range of stellar masses using 
a Salpeter initial mass function (IMF) with a differential 



number of stars, dTV^/dm* oc 



in the mass range 
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O.IMq < m* < IOOMq. The number of disruptions is domi- 
nated by low mass stars despite their smaller rt; switching to 
a top-heavy IM F, as is somet imes discussed in the context of 
galactic nuclei jBartko et al.i,201(l 'l. would reduce Ai'<(100) 
by a factor of a few. 



4 WET MERGERS 

In gas-rich mergers, the pre-kick profile is likely to be differ- 
ent. On the one hand, rapid loss of angular momentum by 
the binary to dynamical friction on the gas can produce a 
core by denying stars the tim e needed to relax in to a cen- 
tral cusp as described above (|Merritt et al.l [20091 '): but on 
the other hand, in situ star formation could rebuild a nu- 
clear cusp while the binary orbit hardens. The possibility 
of star formation, and subsequent migration, in disks moti- 
vates us to consider values of 7 = 1.5, 1.75, 2. Alternatively, 
stars formed el sewhere could be "ground down" into orbits 
inside the disk (|Sver et al.lll99lh . and then behave in a sim- 
ilar fashion. Although the details of star formation in disks 
fragmenting due to gravitational instability are quite com- 
plex ijShlo sman fc Bcgclman 1987; Alexander et al. 2008), 
we provide an approximate description of their potential to 
contribute to the post-kick loss cone here. If the Roche ra- 
dius of the star in the disk exceeds the disk scale height, and 
tidal coupling of the star to the disk is at least as effective as 
viscosity at transporti ng angular mom entum, the star will 
open a gap in the disk (|Sver et al.lll99ll ) and migrate inward 
on a viscous timescale, to the point wher e those conditions 
are n o longer met, or the disk's inner edge (|Goodman fc TanI 
I2OO4! ). whichever is larger. Here we use Equation (O for the 
radius of influence, since we are considering galactic nuclei 
in the process of rebuilding their cusps. 

We follow a similar procedure as with dry mergers, to 
approximate the size of the pre-merger loss cone. Assuming a 
thin disk, the radial size of the central cavity is determined 
by setting t he viscous timescale at the radius of marginal 
self-gravity (jCoodman fc Tanll2004l ). 

rr A r, inS -1/3 -1/2 1/3 / ^0.1 \ j,fl/2 f-,A\ 

Tvis = 4.2 X 10 yr Qq 3' k ' /i ' I — I Mg' , (14) 

equal to Tqw- Here, k is the opacity in units of electron- 
scattering opacity, fi is mean gas particle mass in units of 
the proton mass, ao.3 is the standard (Shakura-Sunyaev) 
viscosity parameter scaled to 0.3, eo.i is the radiative effi- 
ciency scaled to 10%, and Ie is the total radiated luminos- 
ity in units of the Eddington limit for the black hole mass 
Mbh = Ms X 10* M0. Noting the weak power law depen- 
dences in Tvis, we set all parameters except Mbh to their 
fiducial values l|Goodman &: Tanl2004l ). We remove any stars 
with pericenters interior to the radius at which the SMBH 
binary decouples from the accretion fiow, and assume that 
the remaining stars obtain a nearly spherical angular distri- 
bution ( similarly to the in nermost S-stars in the Milky- Way 
nucleus (|Ghez et al ]|2008l )) before the SMBH recoil. Because 
the details of star formation within the disk are highly uncer- 
tain and the TDE rate is dominated by low-mass stars, we 
assume for simplicity ~ IMq in the wet merger case. We 
emphasize that in this case, accretion-induced alignment of 
SMBH spins prior to merger is expected to strongly suppress 
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Figure 2. Expected number of stars disrupted in less than 100 
years, 7V^<(100), in the wet merger scenario. Left panel: Mass de- 
pendence with i)it = 200 km s~i. The solid green Une represents 
a 7 = 1.5 cusp, the dotted blue 7 = 1.75, and the dashed black 
7 = 2.0. Right panel: Velocity dependence with Mbh = IO'^'^Mq, 
and the same lines as in the left panel. 

kicks over 200km s"'^ (jBoedanovic et al Dotti et al.l 

[ioiO). 



5 OTHER CONSIDERATIONS 

Other processes could also partially refill the binary loss 
cone. In analogy to th e problem of resonant capture during 
planetary migration (|Yu fc Tremainel l200ll ). mean- motion 
resonances could be capable of pulling stars inward during 
the final stages of the SMBH merger, dramatically increasing 
the number of post-kick disruptions. The special case of the 
1:1 Lagrange point resonance has recently been investigated 
(|Schnittmanll2O10l ) and found capable of migrating stars to 
wit hin tens of Schwarzs child radii from the system barycen- 
ter (ISeto fc Mutdl2010l ). Higher integer ratio mean-motion 
resona nces have been seen to affect stars (see IChen et al.l 
I2OO9I), Figures 5, 6, 7) as a binary SMBH hardens, although 
their ability to drive resonant migration is less clear. A de- 
tailed study of resonant migration in SMBH binaries is be- 
yond the scope of this paper, but this effect has the potential 
to dramatically expand the short-period population of the 
post-kick loss cone. 

For a source at a redshift 2, cosmological time dilation 
will stretch the duration of each TDE flare, delay the onset 
of the first post-kick flare, and reduce the observed TDE 
rate all by the same factor of {1 + z). However, observations 
suggest that mean density of stars in high-redshift galaxies 
scales as (1 + zf (|Oesch et al.ll2010l ). This could lead to a 
net enhancement in the observed TDE rate cx (1 -I- z)^ per 
galaxy, if the central regions of galaxies are self-similar. We 
set 2 = to ignore these possible cosmological effects in our 
calculated TDE rates. 
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Figure 3. Left panel: Expected number of stars disrupted in 
less than t years, iV<(t), for the dry merger model and A/bh = 
10^-^ Mq. The green hne is for iij; = 200 km s~^, the blue line 
is for Vj^ = 400 km s~^, the red hne is for vy^ = 800 km s~^, 
and the black line is for = 1600 km s~^. Right panel: Af<(t) 
for V]g^ = 400 km and varying masses. The green line is 

for Mbh = 10® M0, the blue for Mbh = lO'^'^M©, the red for 
Mbh = W^Mq, and the black for Mbh = W^'^Mq. 

6 RESULTS 

For simplicity, our calculations assume binaries of equal 
mass black holes. An unequal mass would increase Tqw and 
increase N^{t) by allowing more time for refilling the bi- 
nary loss cone, but it would also decrease the likely val- 
ues of life. The latter change dominates only for mass ratios 
smaller than ~ 0.1, so our results should be regarded as 
conservative for major mergers. Figure [1] presents the veloc- 
ity and mass dependences of our most realistic model (the 
joint core-cusp profile, for SMBH binaries that harden in 
dry mergers by scattering of stars), as well as a less realis- 
tic core-cusp model included for illustrative purposes. The 
first model (labeled BLC) removes stars from the binary 
loss cone and results in an interesting number of TDEs for 
Mbh < lO^M© and 200 km s"^ < Wk < 1000 km s'V At 
low velocities, overlap with the pre-merger loss cone sharply 
suppresses A'^<(100), while at high velocities, the reduced 
size of the bound stellar population also shrinks A^<(100). 
At higher masses, all short-period stars are scoured by the 
pre-merger binary loss cone. The second case in Figure [1] 
(labeled NBLC) replaces the loss cone of a binary with that 
of a single black hole. Interestingly, the dramatic increase in 
A''< (100) here results from the addition of relatively few stars 
(~ 100 for AIbh ~ W^Mq), indicating that resonant migra- 
tion of a small population of stars could significantly boost 
A'^<ioo. Finally, the 7 = 1 (pure core) case does not have 
enough centrally located stars to produce post-kick tidal dis- 
ruptions on < 100 yr timescales. 

The results for wet mergers are illustrated in Figure 
[2] Although these rely on significantly more uncertain as- 



sumptions than the dry merger model, they show that star 
formation can produce interesting values of 7V<ioo in the 
7 = 1.75 and 7 = 2 cases. 

Figure shows how the delay until the first post-kick 
disruptions changes with black hole mass and kick velocity. 
Here we consider the joint core-cusp model (with removal 
of pre-merger loss cone by the binary), and find that the 
first disruption is expected to occur between two and five 
decades after SMBH coalescence for black holes with masses 
between 10^ Mq and 1O^M0, and kick velocities between 
400 and 800 km s~^. Fortunately, this region of parameter 
space falls within both the black hole mass range LISA is 
likely to observe, and the range of physically plausible recoil 
velocities for dry mergers. 



7 SUMMARY 

We find that merging galaxies with black hole masses 
Mbh ^ IO^Mq are likely to produce a tidal disruption flare 
a few decades after coalescence in the case of a gas-poor 
merger. Multiple flares could possibly be seen on a timescale 
of years if resonant migration is effective. For gas-poor merg- 
ers, the peak rate is therefore at least ^ 10"^ times highe r 
than the typical TDE rate in galaxies (|Donlev et al.|[20o3 ). 
The total number of TDEs is maximized, and delay until 
the flrst TDE minimized, if the kick velocity is in the range 
of 2QQ-lQQQkm s'^ 

Our minimal predictions concerning dry mergers could 
be dramatically enhanced if resonant migration increases the 
number and frequency of post-kick tidal disruptions. The 
case of wet mergers is substantially more complicated, and 
because star formation and disk migration can also signifl- 
cantly increase the population of the post-kick loss cone, a 
much more detailed study of star-disk interactions around 
binary SMBHs is needed to make robust predictions in the 
gas-rich case. 

Moderate to high kick velocities in dry mergers will pro- 
vide a robust EM counterpart to the GW signature of black 
hole coalescence within the LISA band, enabling accurate 
identification of the host gal axy and a precise mea surement 
of cosmological parameters ||Ho1z fc Hughes! |2005| ) within a 
few decades of the initial GW signal. With the advent of 
massive transient surveys, such as PTF|f|, Pan-STARRS 0, 
and LSST 0, it is possible that sequential tidal disruption 
fiares could fiag black hole recoil events without a GW sig- 
nal, providing an independent test of the strong field regime 
of general relativity and a calibration of the expected event 
rate for LISA and PTAs. 
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